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Abstract 

The mechanism of hadronic interactions at very high energies is still unclear. Available accel- 
erator data constrain weakly the forward rapidity region which determines the development of 
atmospheric showers. This ignorance is one of the main sources of uncertainty in the determi- 
nation of the energy and composition of the primary in hadron-induced atmospheric showers. In 
this paper we examine the effect on the shower development of two kinds of collective effects in 
high-energy hadronic interactions which modify the production of secondary particles. The first 
mechanism, modeled as string fusion, affects strongly the central rapidity region but only slightly 
the forward region and is shown to have very little effect on the shower development. The second 
mechanism implies a very strong stopping; it affects modestly the profile of shower maximum but 
broadens considerably the number distribution of muons at ground. For the latter mechanism, the 
development of air showers is faster mimicking a heavier projectile. On the other hand, the number 
of muons at ground is lowered, resembling a shower generated by a lighter primary. 

PACS numbers: 13.85.Tp, 13.85.-t, 96.40.De 
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I. INTRODUCTION 



The cosmic ray energy spectrum extends up to several 10 20 eV. The presence of these 
high energy particles (UHECRs) poses an enigma, since many good arguments suggest that 
they should not be observed because of the GZK cutoff This apparent contradiction has 
stimulated a variety of more exotic explanations of their origin and composition which is 
usually known as the GZK puzzle. Both theoretically and experimentally the problem has 
been the subject of considerable efforts [2j. 

One of the major uncertainties is the lack of information about the composition of the 
highest energy cosmic rays. Besides being a difficult, indirect measurement with contradic- 
tory experimental results, it involves extrapolations of hadronic physics models far beyond 
the actual measurements. The systematic effects due to the uncertainty of the hadronic 
interactions, see e.g. (^J, at these energies are far from being established. 

In Cosmic Ray Physics, it is usually assumed that a nucleus-nucleus collision occurs as 
an independent combination of nucleon-nucleon interactions. This approximation is useful 
in that it allows to obtain scaling relations for shower properties like shower maximum and 
the muon content, which are roughly obeyed. But it is not justified at high energies: a beam 
nucleon may interact several times with different target nucleons, which will invalidate the 
approximation trivially. A more interesting possibility would be the existence of a higher 
degree of collectivity (defined in the sense of a failure in the independence of the superposition 
of nucleon-nucleon interactions). This collectivity should be present in any type of hadronic 
collisions, but its enhancement due to the hadron size makes it more prominent for nucleus- 
induced reactions. It could eventually lead to the formation of new nuclear phases during 
the collision, like the Quark-Gluon Plasma. The existence or not of such new form of matter 
at the SPS at CERN and RHIC at BNL has been widely discussed Q]. For example, it is of 
relevance for our study that some violation of the linear scaling law for the rapidity shift of 
the projectile nucleons, which works at lower energies, has already been observed at RHIC 

a. 

Here we will investigate the effects of collective behavior on the development of high en- 
ergy cosmic ray showers. Collective behavior may affect shower development in two different 
ways. First, it could lead to changes in the total cross sections. Such changes have been 
investigated in some models and found to be negligible, see e.g. Second, it is generally 
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expected that any collective behavior would affect the spectra of secondary particles; this is 
the point that we will try to address in our study, exploring it in two directions. On the one 
hand, it may reduce the total multiplicity due to shadowing. As we will see below this would 
affect the shower development only if the multiplicity reduction can be carried over the whole 



cascade process 



,7 



]. This is, in general, not expected. Alternatively, collective effects may 
affect the secondary particle spectrum, in particular in the large Feynman-x region, thereby 
modifying the inelasticity of the collision. This may have important consequences in shower 
development as we will show below. Other possible consequences of collectivity like heavy 
flavor production or modifications of large-pt particle production are expected to have much 
smaller effect on the shower development. 

We will analyze the seco nd p ossibility within both directions. We will consider first the 
possibility of string fusion as a collective effect which strongly reduces the multiplicities 

in the small Feynman-x region but modifies only slightly the large Feynman-x one. It turns 
out that this effect produces negligible modifications in shower development, as expected. We 
then consider a second possibility in which extremely strong collective effects do modify the 
stopping power and inelasticity in the collision. Concretely, it has been suggested recently 
by Mishustin and Kapusta that in the central region of a high-energy collision, strong 
chromo-electric fields are formed. They will attract the forward and backward particles, 
reducing considerably the rapidity of the fast-moving, leading secondaries. This effect is 
expected to increase the inelasticity even for peripheral collisions. It will lead to a small 
change in the shower depth profile but it will induce a large broadening of the number 
distribution of muons at ground. 

The plan of the paper is as follows. In the next Section we will briefly discuss the ingre- 
dients present in hadronic models introduced in the simulation of atmospheric showers. In 
Section III we will analyze the influence of the string fusion mechanism on shower develop- 
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will be discussed in Section IV, and Section V 



ment. The extreme stopping scenario of 
will be devoted to our conclusions. 

II. HADRONIC MODELS 

Most of the models which simulate interactions of nucleons are based on the Gribov-Regge 
Theory (GRT) 12], which most successfully describes elastic scattering and, via the optical 
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theorem, the total hadronic cross section as a function of energy. In GRT the observed rise 
of the cross section at high energies is a consequence of the exchange of multiple supercritical 
Pomerons, the Pomeron being here the effective exchange. Inelastic processes are described 
by cut Pomerons. In the Dual Parton Model (DPM) or Quark-Gluon String Model (QGSM) 
13| . each cut Pomeron leads to two color strings which fragment subsequently into color- 
neutral hadrons. The number of exchanged strings (or Pomerons) rises with energy. Most 
models widely used in Cosmic Ray Physics to compute hadronic collisions at high energies, 
e.g. , 3, 3 , are based on these ideas for the low transverse momentum region; for hard 
particle production, perturbative QCD is used with a matching between the soft and hard 
components. 



A possible formulation of GRT can be done [17( assuming that strings have a transverse 
size A = nr 2 . As the number of strings increases with collision energy, mass number of 
projectile and target and decreasing impact parameter (increasing centrality) of the colli- 
sion, it is expected that the strings overlap and can no longer be considered to evolve and 
fragment independently. The simplest possibility to consider is fusion of strings. The fu- 
sion probability is determined by some parton or string transverse dimension and by the 
corresponding parton-parton cross section. This model is the string fusion model (SFM) 
embedded in the DPM. In the last version of this model, in the form of the program PSM 
[lol | , only fusion of pairs of strings is taken into account and consecutive fusions of strings are 
neglected. Therefore, the string fusion parameter, the transverse area occupied by a string, 
is considered as an effective parameter that is taken to be [10] a v = 2ixr 2 = 7.5 mb. Let 
us note that for proton-nucleus collisions the number of strings is much reduced compared 
to nucleus-nucleus collisions. Therefore even at the highest cosmic ray energies the effect of 
string fusion on proton-air collisions is very small. 

Consequences of string fusion 0, [n]] are the reduction of multiplicities at central ra- 
pidities, a slight increase in the transverse momentum of the produced particles, enhanced 
production of strangeness and of baryon/antibaryons and enhanced correlations between 
particles produced in the forward and backward rapidity hemispheres. Nevertheless, these 
effects are mainly noticeable at central rapidities, the very forward and very backward re- 
gions being only slightly affected [18J. Results of the program PSM for cross sections, 
multiplicities, longitudinal and transverse momentum spectra and production of different 
hadron species can be found in [ijj], together with an extensive comparison to accelerator 



4 



Fe-Air 1 9 GeV 




FIG. 1: Secondary particle energy distribution for Iron-Nitrogen collisions at 10 9 GeV laboratory 
energy without fusion (continuous line) and with fusion (dotted line). 

data. In Fig. ^the energy distribution for secondary particles in Iron- Nitrogen collisions at 
10 9 GeV laboratory energy is shown with and without fusion. The reduction of multiplicities 
for low-energy secondaries due to fusion is clearly seen, while the distribution of high-energy 
secondaries is very modestly affected. 

An appealing possibility to extend the mechanism of string fusion, allowing for a truly 
collective phase, is to consider a percolation phase transition [19] : When the density of 
strings in the transverse space of the collision reaches a critical value, paths of overlapping 
strings crossing the total available area appear and the average cluster size increases suddenly 
with increasing string density. Consequences of percolation have been studied extensively 
j^ . This phenomenon implies the existence of very strong chromo-electric fields stretched 
between projectile and target color sources, like in the model of [uj which will be discussed 
in Section IV. 

As a last comment in this Section, let us note that fusion and percolation of strings 
can be viewed as effective realizations of a collective mechanism of interaction between 
the exchanged Pomerons, as those introduced in 15|, [l6|. This mechanism also leads to a 
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reduction of multiplicities. 



III. AIR SHOWERS AND STRING FUSION 



We have implemented the PSM with/without string fusion (program PSM-1.0 |l0|) in the 



Air shower simulation program AIRES [2l| . In our implementation, the possibility of string 
fusion is considered only for the first nucleus-nucleus collision. This is a valid approximation 
since the energy of the secondary collisions is greatly reduced compared to the initial energy 
and also because these secondary collisions take place mainly between individual nucleons 
and mesons and air nuclei; thus, the density of strings is reduced in the secondary collisions 
making the probability of fusion negligible. The AIRES hadronic model only takes into 
account kaons, nucleons, and pions. All the other hadronic particles produced in the PSM 
(or QGSJET are forced to decay before being passed to the AIRES subroutines. 

We have run our code for iron initiated showers with energies from 10 14 to 10 20 eV and 
with a relative thinning energy (see [2]]]) of 10~ 6 . 500 showers for each energy have been 
generated both with and without string fusion. 

Our results can be seen in Fig. I2HH1 In Fig. |2]we show the A max distribution for showers 
of 10 20 eV, with and without fusion. The difference in the average value of the A max is less 
than 0.6%, compatible with a null result, despite the fact that at these energies the difference 
in multiplicity of the first collision is about a 40%. The average shower development are 
indistinguishable, as can be seen in Fig. |2]for the longitudinal development of electrons and 
in Fig. El for that of the muons. 

These results seem to contradict previous estimations of the effect j?l 0| and also the 
elongation rate theorem [23. However both the elongation rate theorem and the calcula- 
tions of deuces S B -led that the .eduction of multiplicity can be ca„,ed over the 
whole cascade process. For processes like the one discussed here in which the reduction of 
multiplicity is related to the density of interacting strings, this is not expected. From this 
result, we conclude that string fusion (or any collective effect which affects only the central 
rapidity region) has no or very small effect on the shower development. 
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FIG. 2: Average e + e~ longitudinal development for iron showers of 10 20 eV in the case with fusion 
(dashed line) and without fusion (continuous line); also results from QGSJET are shown (dotted 
line). 

IV. AIR SHOWERS AND REDUCTION OF PARTICLE PRODUCTION IN THE 
FRAGMENTATION REGION 

A different possibility that produces a large stopping as a consequence of strong collective 



effects, was recently considered by Mishustin and Kapusta In the center-of-mass frame 
a nuclear collision can be viewed as two thin, Lorentz-contracted, slabs of nucleons colliding 
and going one through each other. The collision proceeds through the exchange of color 
between nucleons of the projectile and target. After the collision the two slabs of nucleons will 
move fast apart, and in the central region a strong coherent chromo-electric field, modeled as 
color strings stretched between the excited color charges in the slabs, may form. The back- 
reaction of this field on the slabs will decelerate them, causing a large stopping power. The 
situation resembles the movement of charged capacitor plates on the electric field created 
by themselves, with a conversion of kinetic into potential energy. This deceleration will 
continue until the field lines stretch so much that they decay into quark- ant iquark pairs 
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FIG. 3: Average p + p longitudinal development for iron showers of 10 20 eV, with the same line 
convention as in Fig. |2I 

and gluons produced from the vacuum via the Schwinger mechanism, which will eventually 
neutralize the chromo-electric field. 

The resulting (projectile or target) slab will have a final 7-factor given by ^j] 



7* = coshy* 



7o 



1 



2A 



(1) 



where 70 = coshyo, ^0 = tanh?/ are respectively the original 7-factor and velocity of the 
slab and tq is the time from the start of the collision until the end of the deceleration process. 
A is the characteristic deceleration length 



A 



£Po 



(2) 



where e is the initial energy per baryon, ej is the energy density of the chromo-electric field, 
/ is the slab thickness depending on the position in transverse plane (impact parameter), 
and po is the nuclear density (p^l being the nuclear profile function normalized to the mass 
number A). Following Reference we parameterize the energy density on the chromo- 
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electric field by 



s\ a/2 /N n N t 



where a ~ 0.3 is given by the low-x structure function behavior reflecting the expected 
increase in the number of produced partons, and eo sets the energy scale, s is the total 
center-of-mass energy squared and sq — 1 GeV 2 is a reference energy squared. N P)t is the 
participant nucleon density of projectile and target, depending on the impact parameter of 
the collision, and iVo — 0.4 fm -2 . The amount of stopping power depends on the initial 
energy and also on the impact parameter. It is expected that €f will be proportional to the 
number of binary parton collisions, which here is parameterized through the dependence on 
the product N p N t . For uncorrelated collisions ~ 1, while in the case of strong correlations 
(like those assumed in percolation of strings |2£|) one expects (3 ~ 0.5 -r- 0.7. 

In this model the results depend on the product €qTq. With these parameters chosen such 
that (e/)ro ~ 6 GeV/fm 2 for central Au Au collisions at center-of-mass energy 65 GeV per 
nucleon, one gets a final center-of-mass energy of 3.5 GeV per baryon [llj]. This can explain 
the strong stopping power already observed at RHIC j^. 

The participant nucleon density N p , N t depends on the impact parameter of the collision. 
The overlapping area of the collision is calculated integrating the profile function of the col- 
liding nucleus. We parameterize the nuclear density using a 3 parameter Fermi distribution 
taken from 2^|. 

The reduction in rapidity can be seen in Fig. H] where we show the 7-factor of the 
secondary nucleons (cosh?/*) as a function of impact parameter for Iron-Nitrogen collisions 
at energy y/s = 200 AGeV. As can be seen in the figure, for small impact parameters the 
reduction factor is large. For central collisions we obtain that nucleons loose more than 60 
% of their energy at = 10 17 eV and more than 80 % for 10 19 eV. 

The effect can be further seen in Fig. El where we show the same 7-factor as a function of 
the energy for fixed impact parameter. This energy loss is compatible with the measurements 
for Au Au collisions made at BRAHMS Q]. 

In Fig. IHlwe show the center-of-mass xf distribution for Fe-Air collisions at 10 18 eV in 
the case of increase and with no increase of the stopping power. The forward and backward 
peaks are clearly reduced. 



We have implemented this increasing of the stopping power in AIRES 



2l| . In a given 
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FIG. 4: Forward secondary nucleon 7-factor (coshy*) as a function of the impact parameter for 
Fe-Air collisions at center-of-mass energy ^/s = 200 GeV per nucleon for the case [3 = 1 (squares) 
and (3 = 0.5 (asterisks). 

nucleus- Air interactions, the impact parameter is chosen randomly and the stopping power 
is calculated according to Eq. ^ The collision is generated with a conventional Monte Carlo 
generator (QGSJET and Sibyll are implemented in AIRES). The energy (and momentum) 
of the leading particles obtained in the collision is reduced to match the stopping power 
predicted and this energy, taken from the leading particles, is redistributed between the 
non-leading particles. Two different version have been essayed, which probably bracket the 
possible range of these models: in version Kl we have made the reduction for all nucleons; 
in version K2, the reduction has been performed on only the participant nucleons of each 
nuclei. For the rest of the showering process the collisions are calculated normally, so that 
only the first nucleus-nucleus collision is modified. For very high energy, E > 10 17 eV, and 
small impact parameters, b < 6 fm, the inelasticity increase due to this mechanism is so high 
that the first collision is indistinguishable from a proton-air collision in terms of inelasticity. 
This has the effect of modifying slightly the average A" max of the shower and the total number 
of muons. 
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FIG. 5: Forward secondary nucleoli 7-factor (coshy*) as a function of the center-of-mass collision 
energy per nucleon for Fe-Air collisions at fixed impact parameter (in fermi as marked) and (3 = 1 
(upper curves) and j3 = 0.5 (lower curves). 

In Fig. [7|we show the distribution of shower maximum for 10 18 eV iron showers. Although 
the average value of X max does not change much, (~ 13 g/cm 2 for the Kl case and ~ 8 
g/cm 2 for the K2 case) the distribution broadens, almost a factor 1.5 for the Kl case. This 
is mainly due to the dependence of the strength of the elasticity suppression on the impact 
parameter. To illustrate this point we show in Fig. [H] the correlation between the average 
X max and the impact parameter of the collision. Collisions with a small first-interaction 
impact parameter produce showers with a lower A" max . Due to the high degree of inelasticity 
of these collisions, these showers develop faster. The effect is further enhanced in the model 
of increased stopping power in which the created fields are large and the stopping power 
is also large. At large impact parameter the effect of this mechanism is negligible and the 
result with and without the modified stopping power is the same. In this way the correlation 
between the first-interaction impact parameter and AT max is increased as can be seen in the 
figure. 

Finally let us turn to muon distributions. In Fig. Elwe show the muon number distribution 
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FIG. 6: Center-of-mass Feynman x distribution of secondary particles in Fe-Air collisions at 10 
eV for the standard QGSJET model (continuous line), Kl (dashed line) and K2 (dotted line). 

at ground. Here a dramatic effect is observed for showers suffering from this effect or not, the 
number of muons resembling that of a proton shower or a regular iron shower respectively. 
Both the average value and the RMS of the distribution change appreciably. In the Kl 
scenario the average number of muons is 7.3 x 10 6 compared to 8.6 x 10 6 for QGSJET iron- 
initiated showers. As a comparison, the average value of the number of muons at ground for 
showers initiated by protons is 7.4 x 10 6 . In the K2 scenario the change in the distribution 
of muons at ground is small. 

In Figure we show the average longitudinal development for the muon component in 
the standard QGSJET and Kl and K2 cases. The muon component for the Kl and K2 cases 
is retarded with respect to the standard QGSJET. 

So we observe two opposite effects of the increased stopping scenario. On the one hand, 
the shower development is faster than in conventional scenarios, which makes the shower 
look as initiated by a heavier primary (i.e. more Fe-like). On the other hand, the number 
of muons at ground is smaller, which makes the shower look as initiated by a lighter pro- 
jectile (i.e. more proton-like). This distorts the N fl -X max relation used conventionally to 
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FIG. 7: Shower maximum distribution for Iron showers at E = 10 18 eV for standard QGSJET 
(continuous line), Kl (dashed line) and K2 (dotted line). 

discriminate composition in an (possibly) energy dependent way. Also it introduces a source 
of systematics in the energy determination based on the muon number such as inclined 
showers. The effect may be as large as the proton-Iron difference in the number of muons. 

V. CONCLUSIONS 

The mechanism of hadronic interactions at very high energies is still uncertain. Existing 
experimental data refer to the central rapidity region, and the forward region is very weakly 
constrained. This uncertainty constitutes one of the main sources of systematic error in the 
study of air showers produced by UHECRs, thus limiting the energy and composition deter- 
mination crucial to clarify their origin. In this paper we have tried to study the consequences 
on shower development of the existence of collective mechanisms in nucleus-nucleus colli- 
sions which modify secondary particle production. These mechanisms have been proposed 
0, Q, 0] to explain several experimental facts in heavy ion collisions. We have tried 
two extreme scenarios which, in our view, contain the most plausible extreme possibilities: 
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FIG. 8: Correlation between first-interaction impact parameter b and shower maximum X max for 
10 18 eV iron initiated showers using the QGSJET model (continuous line) and the stopping power 
enhancement mechanism in the Kl case (dashed line) and K2 (dotted line). 
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first, string fusion [9] which strongly affects multiplicities but has very little effect on the 



forward region; second, a different picture of strong fields 
stopping. More concretely: 



ll| which strongly increases the 



We have implement the program PSM-1.0 which include collective effects like string 
fusion in a Monte Carlo that simulate air showers, AIRES. In this way we have shown 
that string fusion has almost no effect on air shower development because it only 
changes the multiplicity in the central region, the development of air showers being 
determined mainly by the particles in the forward region. 

We have shown that other collective effects like the increasing of the stopping power 
in the fragmentation region can have consequences in the shower development. The 
average X max is moderately reduced; this faster shower development makes the shower 
look as initiated by a heavier primary. The distribution of muons at ground is reduced, 
thus mimicking a lighter primary in conventional scenarios. 
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FIG. 9: Muon number distribution at ground for 10 18 eV iron initiated showers for standard 
QGSJET (thick continuous line), Kl (thick dashed line) and K2 (dotted line). Also shown the 
muon number distribution for proton initiated showers at the same energy (thin continuous line). 

• Correlations between impact parameter and shower properties may arise in these mod- 
els and constitute a novel aspect to be considered. 

These effects may have important implications in the composition measurements and 
separation of proton-induced from iron-induced showers. They should also be taken into 
account to estimate the systematics of the primary energy determination based on shower 
profiles and number of muons at ground. 



VI. ACKNOWLEDGMENTS 

We thank J. Alvarez-Muniz, R. Engel, and E. Zas for suggestions and discussions. This 
work is supported by Xunta de Galicia (PGIDT00PXI20615PR), by CICYT (AEN99-0589- 
C02-02), and by MCYT (FPA 2001-3837). We also thank the "Centro de Supercomputacion 



15 



1.2e+07 



Z 




1000 1200 



FIG. 10: Average longitudinal development for iron showers of 10 eV in standard QGSJET 
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